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ABSTRACT 

"'/In  this  paper  the  basic  semiconductor  device  equations  modelling  a 
symmetric  one-dimensional  voltage-controlled  diode  are  formulated  as  a 
singularly  perturbed  two  point  boundary  value  problem.  The  perturbation 

(JLyuti.  •  ■= 

parameter  is  the  normed  Debye-length  of  the  device.  -We-  derive  the  zeroth  and 
first  order  terms  of  the  matched  asymptotic  expansion  of  the  solutions,  which 
are  the  sums  of  uniformly  smooth  outer  terms  (reduced  solutions)  and  the 
exponentially  varying  inner  terms  ( layer  solutions ) .  The  main  result  of  the 
paper  is  that,  if  the  perturbation  parameter  is  sufficiently  small  then  there 
exists  a  solution  of  the  semiconductor  device  problem  which  is  approximated 
uniformly  by  the  zeroth  order  term  of  the  expansion,  even  for  large  applied 
voltages.  This  result  shows  the  validity  of  the  asymptotic  expansions  of  the 
solutions  of  the  semiconductor  de”'ce  problem  in  physically  relevant  high- 
injection  conditions.  _ 
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SIGNIFICANCE  AND  EXPLANATION 


It  is  well-known  that  potential  distribution  and  current  flow  in  a 
pn-junction  diode  are  described  by  the  solution  of  a  system  of  ordinary 
differential  equations  subject  to  certain  boundary  conditions.  We  scale  the 
system  appropriately  and  obtain  a  singular  perturbation  problem,  i.e.  certain 
derivatives  of  the  dependent  variables  are  multiplied  by  a  small  parameter 
which  is  identified  as  the  normed  Debye-length  of  the  device. 

The  singular  perturbation  character  of  the  problem  introduces  two 
different  scales  of  variation  of  solutions,  namely  a  fast  one  on  which  the 
solutions  vary  close  to  the  pn-junction  and  a  slow  one  away  from  the  junction. 

We  derive  separate  representations  of  the  solutions  which  hold  inside  and 
outside  the  pn-layer  respectively,  and  we  obtain  asymptotic  expansions  of 
solutions  by  matching  these  local  representations. 

The  main  result  of  thiB  paper  is  that  the  asymptotic  expansion 
'represents1  a  solution  of  the  semiconductor  problem  for  small  X,  i.e.  there 
is  a  solution  which  is  approximated  well  by  the  derived  (finite)  asymptotic 
expansions,  provided  the  singular  perturbation  parameter  is  small.  We  present 
an  estimate  for  the  approximation  error  depending  on  the  applied  voltage  and 
show  that  the  singular  perturbation  approach  'covers'  physically  relevant  bias 
ranges  for  modern,  highly  doped  devices. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  KRC ,  and  not  with  the  authors  of  this  report. 


ASYMPTOTIC  REPRESENTATION  OF  SOLUTIONS  OF 
THE  BASIC  SEMICONDUCTOR  DEVICE  EQUATIONS 

Peter  A*  Markovich*  and  Christian  Schmeiser* 

1 .  INTRODUCTION 

This  paper  is  concerned  with  the  asymptotic  representation  of  solutions  of  the  basic 
semiconductor  device  equations  for  the  case  of  a  simple  model  device,  namely  the  symmetric 
one-dimensional  diode. 

The  physical  situation  we  encounter  is  as  follows.  A  semiconductor  (e.g.  silicon)  is 
doped  with  acceptor  ions  in  the  left  side  (p-sids)  and  with  donor  ions  in  the  right  side 
(n-side)  and  a  bias  v  “  “  vc  is  applied  to  the  Ohmic  contacts  (see  Figure  1). 


Figure  1i  Diode 

For  simplicity  we  aurite  the  following  symawtry  assumptions > 

(i)  The  pn-junction  (that  is  the  boundary  between  the  n  and  the  p-region)  is  in  the 
middle  of  the  device. 

(ii)  The  concentration  of  acceptor  atoms  in  the  p-side  and  the  concentration  of  donor 
atoms  in  the  n-side  are  constant  and  equal  (to  c  >  0,  C  is  called  doping  concentration). 

(iil)  The  applied  potentials  VC,VA  satiafyi  Vc  ■  -VA. 

Under  these  assumptions,  the  performance  of  the  device  is  described  by  the  following  two- 
point  boundary  value  problem 
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A  t  -  n  -  p  -  1 


n*  -  n*1  +  ^ 


»*  .  -nil1  -  i 


0  <  x  <  1 


subject  to 


( 1 • 1e) 


♦  (0)  -  0,  *<  1)  -  *  (6)  -  ♦  (6)  -  tn(  —  /l  *  46  ) 


(l.lf) 


n(0)  -  p(0),  n(  1 )  -  -j-  (l  +  /l  +  4«4J 


(l.lg) 


P(1)  -  4  (-1  ♦  A  +  4«4) 


where  the  dependent  variables  have  the  following  physical  weaning 
♦  t  (scaled)  electrostatic  potential 
ni  (scaled)  electron  concentration 
pi  (scaled)  hole  concentration 
J:  (scaled)  total  current  density. 

The  x-lnterval  [0.1]  on  which  the  problea  is  posed  represents  the  n-side  of  the  device 
v  1 

(after  scaling),  0  »  —  (where  ■  -jj-  Volt  is  the  thermal  voltage)  represents  the 
UT 

voltage  parameter,  is  that  potential  at  x  -  1  which  prevails  if  sero  external  bias 

is  applied  (built-in-potential,  originating  from  the  doping)  and  1,6  >  0. 

The  problea  (1.1)  is  derived  by  adapting  the  basic  seaiconductor  device  equations  (as 
given  by  Van  Roosbroeck  (1950))  to  our  specific  device  using  the  syaaetry  assumptions  (i), 
(ii),  (iii)  and  other  simplifying  assumptions  (like  constant  electron  and  hole  mobilities 
and  neglect  of  generation-recombination  of  carriers)  and  by  appropriate  scaling  (the  length 
of  the  device  is  scaled  to  2,  the  doping  concentration  to  1).  Details  on  the  assumptions 
and  on  the  scaling  are  given  in  Markovich  and  Ringhofer  (1984)  and  Markovich  (1983), 

(1984). 
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The  ptramteri  X  and  6  are  Introduced  by  the  scaling.  Physically  X  is  the 
noraed  Debye  length  of  the  device,  that  means 


where  £  is  the  material  permittivity  constant,  q  the  elementary  charge  and  2t  the 

(original)  length  of  the  device. 

2 

6  is  the  scaled  intrinsic  number  of  the  device 


2  ni 

(1.2b)  *  C 

(nA  is  the  number  of  free  electrons  of  the  semiconductor  per  unit  volume,  n^  “  1010/cm3 
for  silicon). 

e,UT,q  and  nA  are  material  constants  while  the  length  l  and  the  doping 
concentration  C  specify  the  device. 

For  a  realistic  silicon  diode,  we  have  C  >  1017/cm3  and  £  «  10“2cm.  This  yields 
X  <10“'.  Thus  the  problem  (1.1)  can  be  regarded  as  singularly  perturbed  two-point 
boundary  value  problem  with  perturbation  parameter  X  and  the  standard  method  of  matched 
asymptotic  expansion  can  be  employed.  It  turns  out  that  a  boundary  layer  (i.e.  a  small 
region  of  fast  variation)  occurs  in  *>,  n,  p  at  the  junction  x  ■  0.  No  layer  occurs  at 
the  Ohmic  contact  x  *  1,  those  solutions  which  are  approximated  by  the  asymptotic 
expansion  are  uniformly  (in  X)  smooth  away  from  x  «  0. 

Also  S  is  small  (practically  <  10"7).  Me  will,  however,  see  later  on  that  the 

smallness  of  6  has  a  weaker  impact  on  the  solution  structure  than  the  smallness  of  X 
2 

(6  is  not  a  singular  perturbation  parameter  since  it  does  not  multiply  a  derivative  of  a 
dependent  variable). 

In  the  recent  past,  many  papers  dealt  with  the  singular  perturbation  analysis  of  the 
basic  semi  conductor  device  equations  (e.g.  Vasileva  and  Butuzov  (1976),  Vasileva  and 
Stelmakh  (1977),  Markowich  and  Ringhofer  (1984),  Markovich  (1984)).  The  authors  of  these 
papers  concentrate  on  deriving  the  zeroth  order  terms  of  the  asymptotic  expansions  of  the 


solutions  of  the  semiconductor  device  equations  ss  1+0+  (even  in  the  high  dimensional 
esse  and  for  devices  much  more  complicated  than  our  symmetric  diode).  However,  the 
question  of  the  validity  of  these  asymptotic  expansions  is  a  rather  unsettled  issue.  It 
was  proven  by  Markowich  (1984)  that  for  sero  applied  bias  (0  *  0)  there  is  a  solution  of 
the  send. conductor  device  problem  which  is  close  to  the  zeroth  order  term  of  the  expansion 
if  X  is  small  (even  in  the  multidimensional  case).  For  the  one-dimensional  case  this 
result  was  carried  over  to  small  applied  bias  (ju|  <  c  where  c  tends  to  zero  rapidly  as 
X  +  0+,  6  *  O+i  see  Markowich,  Ringhofer,  Selberherr  and  Langer  (1982)).  However,  this 

is  of  extremely  limited  practical  applicability,  since  it  only  means  that  for  sufficiently 
high  doping  or  for  a  sufficiently  large  device  ( large  C  or  large  l  implies  small  X ) 
and  for  biasing  conditions  sufficiently  close  to  thermal  equilibrium  (represented  by 
0*0)  the  solutions  are  asymptotically  represented  by  the  asymptotic  expansions.  In 
practice,  however,  one  is  interested  in  the  performance  of  highly  doped  devices  when  high 
voltages  are  applied. 

In  this  paper  we  show  -  at  least  for  our  simple  model  device  -  that  there  is  a 
solution  of  (1.1)  which  is  approximated  by  the  zeroth  order  terms  of  the  asymptotic 
expansion  if  X  is  sufficiently  small  and  if  U  is  in  some  specified  voltage  range  which 
gets  larger  as  X  +  0+  and  which  includes  physically  relevant  high-injection  biasing 
conditions.  This  result  in  some  sense  justifies  the  singular  perturbation  approach  for  the 
semiconductor  device  equations. 

The  paper  is  organised  as  follows.  Section  2  contains  the  derivation  of  the  zeroth 
and  first  order  terns  of  the  asymptotic  expansions,  the  main  result  of  the  paper  and  the 
functional  analytic  framework  employed  for  its  proof  are  given  in  Section  3.  The  highly 
technical  details  of  the  proof  are  collected  in  the  Appendices  A  and  B. 


We  will  now  apply  the  approach  of  matched  asymptotic  expansions  to  the  problem 

(1.1) .  Therefore  we  assume  that  the  solution  (f,n,p,J)  has  an  asymptotic  expansion  in 
powers  of  1  where  each  coefficient  is  the  sum  of  a  term  which  is  independent  of  X  and 
of  boundary  layer  terms.  It  was  shown  by  Markovich  and  Ringhof er  (1984)  that  the  terms 
representing  the  layer  at  x  -  1  (at  the  Ohmic  contact)  vanish  in  seroth  and  first 
order.  Thus  we  expect  a  layer  only  to  occur  at  the  junction  x  -  0  and  make  the  ansatz: 

*(x,X)  -  *(x)  +  *(t)  +  X$j(x>  +  (x)  +  ... 

n(x,X)  -  n(x)  +  n(T)  ♦  Xn^(x)  +  Xn^T)  ♦  ... 

(2.1)  . 

p(x,X)  -  p(x)  +  p (1 )  +  Xpj(x)  Xp^d)  +  ... 

j(x)  -  3  +  x51  +  ... 

where  the  dots  denote  terms  of  order  at  least  X2.  Note  that  the  current  density  J  is 
independent  of  x  because  of  (l.ld).  T  “  y  is  the  fast  independent  variable.  The  terms 
marked  with  ' are  independent  of  X  and  the  terms  marked  with  '* 1  are  zeroth  and 
first  order  layer  terms  which  are  required  to  decay  to  zero  as  the  fast  independent 
variable  tends  to  infinity i 

(2.2)  ♦(*)  "  i^**)  *  n(“)  *  n^m)  “  p(“)  *  pjl*)  “  0  . 

By  inserting  (2.1)  into  (1.1)  and  by  equating  coefficients  of  equal  powers  of  X  we  will 
obtain  boundary  value  problems  for  the  terms  in  (2.1). 

Construction  of  zeroth  order  terms i 

By  comparing  coefficients  of  X°  and  evaluating  away  from  x  »  0  we  obtain  the 
reduced  equations 


0  »  n  -  p  -  1  , 


or  equivalently 


n  ■  p  +  1 


2p  *  1 


2(2p  ♦  1> 


For  the  zeroth  order  boundary  layer  terns  we  get  (by  comparing  coefficients  of  X  close 


to  x  ■  0) 


♦  “  n  -  p 


n  -  (n(0)  ♦  n)p  , 


p  “  ~(p(0)  p)$  , 


where  the  dot  denotes  differentiation  with  respect  to  T.  Equating  zeroth  order  terns  in 
the  boundary  conditions  (l.le)  -  (l.lg)  gives 

?(0)  +  *<0)  -  0,  n(0>  +  n(0)  -  p(0)  +  p(0)  , 


*<1>  -  +BI  -  f.  n(1)  -  p,  +  1,  p<  1)  -  p,  . 

AAA 

By  conditions  (2.2),  ♦,  n,  p  have  to  be  trajectories  on  the  stable  manifold  of  (2.4) 

which  is  represented  by 


n  -  n(0)(eT  -  1)  , 


p  -  p(  0)  (e'*’  -  1) 


♦  “  -  *  2(  n  +  p  -  sgn  i(i(0)  , 


(see  Schneiser  and  Weiss  (1984)). 


Using  (2.6)  we  obtain  the  reduced  boundary  conditions  from  (2.5) 

(2.7)  S(0)e"*tQ1  -  ^(0)e^(0),  *<U  -♦„-!»  PM)  -  P,  • 

Equations  (2.3)  can  be  integrated  and  the  solution  of  (2.3),  (2.7)  (reduced  solution)  can 
be  written  as 


•  f.  *  .  •  .  ■ 


*  a  •  _  •  _  »  ,*  ,  e  .  •  m 
%*  »,*  %*  '/  * 
v  ^ **./  .•*./ 

Ovv/.Xvin'v. 


(2.8) 


5  -  £  H  ♦  /(B-f)2  "  25x)  ' 


n  -  p  +  1  , 


$  -  2p  +  in  /b  -  Bt  ' 


where  the  parameter  b  >  1  is  related  to  U  by 

(2.9)  ~  -  ♦gj  +  2Pl  -  -  in  /b  «i  F(b) 


Obviously 

lim  F(b)  -  ",  lim  F(b)  -  F’(b)  - - - - -  -  4r  <  0.  be(1,«) 

b+1  b+"  (b  -  1)Z  2b 


holds.  Thus,  (2.9)  is  a  one-to-one  relation  between  U  e  (-",")  and  b  e  (1,">. 

From  (2.8)  and  (2.9)  we  get  a  first  order  (formal)  asymptotic  approximation  of  the 
jltage-cur rent-character is tics 


£  -  ln  LI. -  471?  ♦  4  ♦  46  4  ♦  2J 
2  4?T1£ 


He  compute  n(0),  p(0)  in  terms  of  b  by  using  (2.8)  and  rewrite  (2.6)  as 

*  b  * 

(2.1ta)  n  -  b  -_-j-  (er  -  1)  , 

*  1  -V 

(2.11b)  P  -  £™-j  (e  -  D  , 


(2.11c) 


b  +  1\ 
b  -  V 


(2. lid) 


*<0) 


$(0) 


-in  /F  . 


(2.11c)  is  obtained  by  observing  that  sgn  <H0)  -  -1  (since  sgn  ?(0)  -  1  and 
f(0)  *  i>(0)  ■  0).  The  unique  solvability  of  (2.11c,d)  and  the  exponential  decay  of  ♦,  n 
and  p  follow  from  an  application  of  the  theory  of  Fife  (1974)  to  the  second  order 
equation  in  (2.4)  (see  Karkowich,  Ringhofer,  Selberherr  and  Langer  (1982)). 


V  rs. 


Construction  off  first  order  coefficients: 

Generally  n-th  order  coefficients  of  Hatched  asymptotic  expansions  are  defined  by 
problems  which  are  linearized  versions  of  the  problem  defining  the  zeroth  order 
coefficients  with  inhomogeneities  only  depending  on  coefficients  of  order  at  most 
(n  -  1)  (see  Schmeiser  and  Weiss  (1984)).  The  first  order  coefficients  satisfy 


",  “  P,  » 


(2.12) 


1  -  .  2J 

J1  +  -  2  P1 
2p  +  1  (2p  +  1) 


2(2p  ♦  1)  <2p  +  ir 


for  the  part  which  is  independent  of  1  and 


*1  '  "I  -  p1  ' 


(2.13) 


b  * 


e  *1  +  (n,  +  tej(O)  +  n** (0)  ♦  VTp’(O)  , 


P,  -  -  e“S,  "  CP,  -  <>5,(0)  -  pj'(0)  -  tTp'(O)  , 


for  the  layer  terms.  The  boundary  conditions 


(2.14) 


3,(0)  +  i,(0)  -  0,  n,(0)  +  0,(0)  -  p,( 0)  ♦  p,(0) 


*,<  1 )  -  0.  11,(1)  -  0,  Pl(1)  -  0 


hold.  The  results  of  Appendix  A(a)  imply  that  decaying  solutions  of  (2.13)  can  be  written 
in  the  form 


■*.  •  **,  *  *  *  >,■  • 
.* 


*v 


V(0) 


(2.15) 


«  — £ —  +  p1(0)(ey  -  1)  +  n 


b  -  1 


r1  b 


i  -  -j, 

-r-f  e  +  p,(0)(e  v  -  D  +  p  , 


where  c  is  a  constant  and 


n  ■»  -p’  (0 )  j  e 


*(x)-*(s) 


(2n(s)  +  s<i(s))ds 


(2.16)  p  -  -p’(0)  /  e 
P 


-H»(T  )>^(S  ) 


( 2p  ( s )  +  si|i(s))ds 


T 

X 


♦  »  — 4>tT )  J  ♦“  (  s)  /  f(tt>[p1(0)(eVl  -  e  V  U')  +  n  (u)  -  p  (u)]duds 


Integration  of  (2.12)  using  (2.14)  and  (2.15)  gives  expressions  for  the  reduced  solution  of 
order  1 


3i  ■  ittt  <V0) '  V0>>  ' 


(2.17) 


1  -  x 


n1  “  P1  “  J1  *1  “  2p1  ’ 


2(2p  +  1) 


The  constant  in  (2.15)  is  c  •  — — r— ^ -  (p  (0)  -  n  (0)).  The  expressions  (2.15),  (2.17) 

(b  +  1 ) 2  P  p 


can  be  simplified  using  the  following 


Lemma  2.1S  n_  ( 0 )  -  p_  ( 0 )  -  -  - —  i(i(0) 

-  p  p  2/b 


Proof i  Relations  (2.3)  and  (2.8)  imply 


p'(0) 


J(b  -  1) 


2(b  +  1) 


Thus,  we  get 


3.  TUB  RtPMBKHTATIOH  THEOREM 

The  question  to  be  answered  now  Is  whether  and  how  well  do  the  zeroth  order  terms  of 
the  asymptotic  expansions  (2.1)  approximate  a  solution  (i|r,  n,  p,  J)  of  the  problem  (1.1). 

The  following  theorem  is  derived  by  a  straight  forward  application  of  the 
representation  theory  for  singularly  perturbed  two-point-boundary  value  problems  given  in 
Schmeiser  and  Weiss  (1984). 

Theorem  3.1a  For  every  U  e  R  and  6  >  0  there  is  Xq  -  X^ (u ,6 )  >  0  such  that  for 
X  <  Xq(U,6)  there  is  a  locally  unique  solution  (^i,  n,  p,  J)  of  (1.1)  which  satisfies 

*(x,X)  -  +  (x)  +  *(£)  +  0(X) 

n(x,X)  ■  n(x)  +  n(j)  +  Q(X) 

p(x,X)  •  p(x)  +  p(^)  +  0(1) 

J  «  3  +  0(X) 

uniformly  for  x  €  (0,1]. 

Also  the  existence  of  asymptotic  approxisiations  of  arbitrary  order  of  accuracy  can  be 
concluded  from  Schmeiser  and  Weiss  (1984). 

The  theorem  says  that,  for  given  doping  concentration  C  (which  determines  6  by 
(1.2b))  and  for  a  given  voltage  U  there  is  1q  >  0  such  that  for  all  diodes  with  length 
21  >  2t g  the  corresponding  semiconductor  device  equations  (1.1)  have  solutions  which  are 
uniformly  approximated  (to  order  X)  by  the  zeroth  order  terms  of  the  asymptotic  expansion 
(1.1)  (note  that  t  and  C  determine  X  by  (1.2a)  for  a  given  material  and  that  X  ♦  0 
for  1  ♦  »  when  C  is  constant). 

Practically,  this  is  a  rather  meaningless  result  since  one  is  not  interested  in 
applying  a  fixed  voltage  to  various  (sufficiently  large)  devices,  but  in  varying  the 
voltage  applied  to  a  fixed  highly  doped  device  (that  means  X  and  5  are  small  and  fixed 
and  U  varies  over  a  certain  voltage  range). 


To  prove  an  approximation  result  which  is  uni fore  in  a  sufficiently  large  0  interval 
we  basically  proceed  as  Schmeiser  and  Weiss  (1984)  did  for  general  singularly  perturbed 
two-point  boundary  value  problems  but  we  always  keep  track  on  how  small  A  has  to  be  (in 
dependence  of  6  and  the  parameter  b  which  is  related  to  U  by  (2.9))  in  order  to 
guarantee  the  validity  of  the  results  in  Schmeiser  and  Weiss  (1964). 

We  regard  the  problem  (1.1)  (after  eliminating  (l.ld),  putting  all  terms  on  the  left 
hand  side  and  expressing  U  as  function  of  b,  $  by  using  (2.9))  as  operator  equation 

(3.1a)  ,6  rb^'l''n'I>' “  <X2*'“  *♦'-  §.  P,+  P*'  +  J'  *<0).*<1)  - 

n(0)  -  p(0),n<1)  -  j  (l  +  /l  +  454),p(1)  -  y  (-1  +  /l  ♦  464JJ  -  0 
(3*1b)  FA,6,b  !  B1,X  B2,A  ‘ 

The  following  spaces  will  be  needed  in  the  sequels 

(3.2a)  B1  A  “  w2'"<°«1)  *  <C10,1]  n  W1,1<0,1))2  x  R 

(3.2b)  B2  A  “  L"(0*1)  x  tl*1  <0,1 )  )2  x  R5 

equipped  with  the  norms  I • I 1  ^ ,  I • I 2  ^  resp.  defined  by 

(3.3a)  l(*,n,p,J)l  .  s-  1*1  „  +  Al*’l  „  +  A2I*"I  „ 

I>  (0,1)  L  (0,1)  L  (0,1) 

+  Ini  „  Ain' I  ,  +  Ipl  „  +  Alp' I  +  |j| 

L  (0,1)  L  (0,1)  L  (0,1)  L  (0,1) 

for  (*,n,p,J)  e  Bj  j  and 

(3.3b)  l(u,v,w,a)l2  ,  «■  lul  m  +  Ivl  .  ♦  Iwl  .  +  la^ 

L  (0,1)  L  (0,1)  L  (0,1) 

for  (u,v,w,a)  e  B2 

WM*P(0,1)  denotes  the  Sobolev  space  of  real  valued  functions  defined  on  (0,1)  whose 
weak  derivatives  of  order  up  to  m  are  p-integrable  for  1  <  p  <  •  and  essentially 
bounded  for  p  -  «  (see  Adams  (1975) )»  i.e.  f  e  tf“'P<0,1)  iff  f(i)  e  lP(0,1)  for 
0  <  i  <  a.  We  denoted 


(3.4) 


Ifl  _  «-  (J  |f(x)|pdx)1/p,  1  <  p  <  -l 

I.P(0,1)  0 


If I  m  l»  sup  jf (x) I  . 

l  (0,1)  xeto.i) 


and  |o|  i-  MX  |a  |  for  a  ■  (a.,..., a  )  e  Rn.  Also  ths  lnducsd  matrix  norm  Is 
KKn  1  n 

dsnoted  by  |*|— . 

F roa  sobolsv's  imbedding  theorem  (see  Adams  (1975))  wo  conclude  that  f  e  W1 ' 1 ( 0, 1 ) 

implies  f  e  CB(0,1)  (f  is  continuous  and  bounded  in  (0,1))  and  that  Ifl  m  < 

L  (0,1) 

const  (ifl  |  +  Ifl  .  J  holds  for  all  f  e  W*'^(0,1).  Therefore  the  functional 

L  (0,1)  L  (0,1) 

Ifl  „  ♦  Xlf'l  .  defines  a  norm  (for  1  >  0)  on  C[0,1]  n  w1r1(0,1)  under 

L  (0,1)  L  (0,1) 

which  this  space  is  a  Banach  space. 

We  denote  the  asymptotic  expansion  (2.1)  (up  to  the  first  order  term)  by 


(3.Sa) 

♦4(x,X) 

(3.5b) 

na(x,X) 

(3.5c) 

P4(*,X) 

( 3 . 5d ) 

J4(X)  - 

-  ¥u)  ♦  ♦  ij^x)  + 

-  ii(x)  ♦  n(y)  *■  Xn j(x)  +  Xn^j^) 

-  p(x)  ♦  p(^)  +  Xpf(x)  ♦  Xp,Cx3 

5  ♦  15,  . 


* 

The  terms  on  the  right  hand  side  of  (3.5>  were  computed  (up  to  ♦  which  solves  (2.11c,d)) 
in  Section  2.  Clearly  n>(  p#  and  J€  depend  on  5  and  on  the  parameter  b. 

For  the  representation  proof  we  proceed  as  followst  (A)  We  compute  the  residual  at 

<*m#Ba'pa'Ja> 

(3.6)  r(X,6,b)  «  F^  j^^(1>a»na,pa'Ja* 

and  estiMte  lr(X,6,b)l_  .  «  R(X,6,b). 

(B)  We  compute  the  Frechet  derivative  of  F,  .  .  at  (<>-,n.,p  ,J  ): 

A  f  O  f  D  A  AAA 


(3.7) 

and  estiMte  its  inverse 


°(t,n,p,J),X,6,b*^a'na'pa'Ja*  *  B1,X  *  B2,X 


KD 


(f.n.p.Jj'x ,«  ,b‘ W'a'V 1  1 B  ^ 


«  K^(X,6  ,b) . 
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(C)  Ns  sstimsts  ths  Lipschlts  constant  of  D.  J.  •  .  i  B«  .  ♦  L(B.  t  ,  )  in  s 
sufficiently  large  sphere  about  (♦^.n^.p^.J^H 


<3*8>  *Dt»n,p,  JI'x,6,b^'n'p'3*  "  D*.n,p.jrX.«.b<VVpa'Ja,,B1  X*B2  , 

<  K2<M,b)l<*  -  ♦*.»  -  n>#p  -  pa,J  -  Ja>«1|X 


for  •  ( v  "  9^/B  "  na,p  -  p#/J  -  *  P* 

Th«  bound*  R,Kj  and  Kj  ara  explicitly  determined  in  the  following 


3.1.  There  ara  Constanta  Xfl  »®o ,C1 ,C2  ,C3 >  0  auch  that  the  ••tiaates 


(3.9) 


R(X,«,b>  <  c.X2l - - - r  ♦  (tn  -fj«8(in  /b)’V2J 

1  (b  -  1) 


Kj(X,«,b)  <  -f- 


(3.10) 

hold  indapandantly  of  the  radiua  p  for  0  <  i  <  Xj,  0  <  t  <  Sj.  If,  additionally 

>  « . .  .  a 

(3.11) 


X(tn  yH - - — nr  +  4Bb<*“  •'b)19/2J  «  a 

(b  -  1)  8 


than 


(3.12)  K  (X,«,b)  <  cJ - - - t+  (in  /F)2J 

(b  -  1)° 

holds,  too. 

Proof i  The  estimate  (3.9)  is  provan  in  Appendix  B,  (B.1)>  (3.12)  in  Appendix  A,  Leans  A.4 
(which  uses  (3.11))  and  (3.10)  follcws  iaaediately  by  linearising  rx  6  b  at 

lVVpa,Ja)- 

He  will  eaploy  the  following  version  of  the  lap licit  function  theorem  (sea  Chow  and 
Hale  (1982)) i 

Lewaa  3.2.  Let  r  be  a  Upechits-contlnuously  Prachat  differentiable  nap  froa  the  Banach 
spaca  B1  into  the  Banach  space  Bj.  Assuan  that  the  aquation  F(x)  -  r  has  a  solution 
xa  which  is  isolated,  i.a.  the  Frechet  derivative  DxF(x4)  t  B,  ♦  Bj  is  boundedly 
invertible.  Let  the  constants  K<  and  Kj  satisfy 


,tD««v>  v*,  *  ki  • 

IDXF(X>  -  Oxr(xa)lB^B  <  Kjlx  -  x#«B^  if  IX  -  x^l  <  p  for  iom  p  >  0 


than  the  problea  F(x)  -  0  ha a  a  solution  x  which  satisfies 


lx.  -  xl.  <  »,lrl_ 
a  B,  1  Ba 


(3.15) 


*B  *  — 5 — ) 


1  2K,Kj 


holds.  Tha  solution  x  is  uniqua  in  tha  aphara  with  radius  ,r  ,  and  cantar  x,. 
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Tha  main  result  of  tha  papar  is  tha  following! 

Thaoraa  3.1.  Thara  «ra  constants  ^ ,C+,C_,D+,D_  >  0  such  that  tha  problaa  (1.1)  has  a 

locally  uniqua  solution  (F.n,p,J)  for  X  <  X^,  S  <  if  either 
«i»  |>Sl  X(|  -  ^  *  l)"  <  C+ 


<ii,  *<♦„  - 

holds.  Tha  solution  (*,n,p,J)  satisfies  tha  aatiMta 


<♦-(♦+  +>l  _  +  In  -  (n  ♦  n)l  m  +  Ip  -  (p  +  p)l  „ 

L  (0,1)  L  (0,1)  I.  (0,1) 


♦  J  -  J  < 


D+X(j  -  tBI  +  l)1/2  in  the  case  (i) 

D  X52(^bi  -  j  +  l^e-0^2  in  the  case  (ii) 


Proof :  A  simple  calculation  which  uses  (2.9)  shows  that  (i)  and  (ii)  imply  (3.11). 
Therefore  the  estiMtas  (3.9),  (3.10)  and  (3.12)  of  Leaaaa  3.1  hold  and  the  implicit 
function  thaoraa  (Leans  3.2)  can  be  applied  to  conclude  the  locally  unique  solvability  of 
Pi  a  w(*,n,p,J)  -  0.  It  follows  that  thara  is  a  locally  uniqua  solution  (f,n,p,J)  if 
r(A ,£,b)  «  — = - - - .  This  is  guaranteed  by  (3.11),  too.  The  inequality  (3.16) 


2Kj(X,4,b)Kj(X,6,b) 


is  obtained  by  aatiMting  tha  first  order  tanas  of  tha  axpansion  (2.1)  (sea  Appendix  A, 


1  .  ■  .  •  Z — '.1  11 — ^ — : — : — ■ — r~.  — ,-T;  _sr  w  m _ 2 _ ! 

*.  *•  *,*  *  '  •  *  ■  *  a  *  •  *  *  ***,  **,  •'  •*  «*  »*  ,*  »•  ‘  •  P  •  ’  a  *  •  ‘  .  *  ,**  '«  ,*.  ‘ •  *,  V  «  *  • 

*.  ,%  •  */  •;  /  %*  V  a*  -  -  .  •  a  -  .  •»•*,  a*  a*.  V  V  V  -  *  »**  •  *  .*•  .*•  ,%  ,*•*  \*  *.*  - 


L»— ■■  A.d.  16,  A.d.  20)  and  by  observing  that  tha  error  bound  given  by  (3.14)  (which 
includes  the  first  order  teres)  is  sharper  than  (3.16). 

He  now  consider  a  realistic  silicon  diode  with  the  numerical  values  of  the  parameters 
given  in  Section  1.  The  cases  a)  and  b)  in  tha  formulation  of  Theorem  3.1  yield 
approximate  upper  and  lower  bounds  for  the  applied  voltage  V  -  Uj>U.  with 
C+  -  C_  «  0(1)  we  obtain 

(3.17)  -0,2  Volt  £  V  £  0,8  Volt 

The  upper  bound  represents  a  large  forward  bias  for  the  considered  device.  Thus  our  theory 
covers  a  realistically  large  forward  bias  range. 

In  the  reverse  bias  case  application  of  up  to  -10  Volt  is  of  practical  interest. 
Apparently  the  reverse  bias  range  covered  by  the  presented  theory  is  such  smaller  than 
desired,  we  conjecture  that  this  is  caused  by  the  fact  that  tha  space  charge  region 
wideness  with  increasing  reverse  bias.  This  limits  the  validity  of  an  ansatz  which 
strongly  uses  the  'boundary  layer  behavior*  of  the  solution. 

Figure  3.1  shows  the  reduced  potential  i.  Figure  3.2  the  reduced  hole  concentration 
p  and  Figure  3.3  the  reduced  current  density  3  as  given  by  (2.0).  The  surfsces  were 

obtained  by  parametrising  the  curves  if .  n  on  the  x-interval  [0,1]  (n-side  of  the 

u  u 

device)  with  the  scaled  applied  voltage  —  ranging  from  “  ■  -10  (-0.5  volt  reverse  bias) 

to  y  «  40  (2  volts  forward  bias).  The  reduced  electron  concentration  n  is  not  depicted 
since  it  is  given  by  shifting  p  by  the  value  1,  i.e.  n  -  p  +  1.  The  extension  of  the 
reduced  solutions  to  the  interval  [-1,1]  is  obtained  by  using  the  symmetry  conditions 

n(x)  -  p( —x ) ,  0(x)  “  -5(-x)  . 

The  reduced  current-voltege  characteristic  3  “  J(U)  exhibits  the  well-known  exponential 
behaviour  (see  Sse  (1981)). 

The  lay er-eguat ions  (2.11)  were  solved  numerically.  Figures  3.4,  3.5,  3.6  show  the 

*  *  *  o 

layer-terms  ♦  <*)»  n(t)  and  p(i)  reap.,  again  as  surfaces  parametrized  by  —  varying  in 

the  range  specified  above.  Figure  3.4  demonstrates  the  increase  of  the  layer  jump  and  of 

the  width  in  the  potential  9  «•  the  reverse  bias  increases.  Figures  3.5,  3.6  show  the 
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depletion  of  carriers  occuring  in  the  space-charge-region  for  reverse  bias  and  the  increase 
of  the  carrier  concentration  in  the  space  charge  region  for  forward  bias. 

Also,  the  'full'  singularly  perturbed  problem  (1.1)  was  solved  numerically  (by  using 
the  general  purpose  two-point-boundary  value  problem  COLSYS  authored  by  Ascher, 

Christiansen  and  Russell  (1976)).  Figures  3.7,  3.8,  3.9,  3.10  depict  the  potential  p, 
the  electron  concentration  n,  the  hole  concentration  p  and  the  electric  field  V' 


A  comparison  of  the  full  solutions  p,n,p  and  the  reduced  solutions  P,n,p  clearly 
shows  that  p,n,p  are  approximated  well  by  P»n,p  outside  the  layer-region  on  the  whole 
considered  bias-range,  in  fact  the  corresponding  reduced  and  full  solutions  agree  at  least 
up  to  plot  accuracy  outside  the  layer.  The  current-voltage  characteristic  J  •  J(U)  is 
not  depicted  since  it  is  graphically  indistinguishable  from  the  reduced  current  voltage 
characteristic  J  “  J(U)  as  shown  in  Figure  3.3.  Note  that  the  depicted  bias  range  is 
larger  than  (3.17). 

Figure  3.10  demonstrates  the  occurrence  of  a  very  large  electric  field  within  the 
layer  (note  the  scale  on  the  PSIP-axisI),  which  is  explained  by  the  asymptotic  expansion 


for  p ' : 


P'(x,X)  ~  V(x>  +  x  + 


■  V  •  •  •  '  .  *\ 


’  >  ’  *  ’  i,'*  ,**  ** 


.*  w  *\  '  . 
,  ■*  * » 

*  - 

•  A  .' 


a 


ELECTRIC  FIELD 
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APPENDIX  A:  THE  LINEARIZED  PROBLEM 

We  shall  analyze  the  linearization  of  problem  (1.1)  at  the  foraal  approximation  of  the 
solution,  which  has  been  constructed  in  Section  2.  Thus  we  consider 

A2?"  -  n  -  p  ♦  at  , 

a  e 

n*  -  (*•  +  y  *  +  A${  ♦  ♦1)n  +  (n  ♦  n  +  ln1  *  An^**  +  J  +  «2  » 

(A.1)  . 

p*  .  +  |  (i  +  A^*  +  *1]p  -  (p  ♦  p  +  XPl  ♦  Ap^  )t*  -  j  +  a3  , 

*(0)  -  a4>  n(0 )  -  p(0)  -  a5,  ♦( 1 »  -  a6,  p(1)  -  a?,  n(1)  -  a0  , 

with  y  *  (<|>,n,p,J)  e  ^  and  a  -  («■),. ..,a0)  e  Bj  We  proceed  as  in  Schmeiser 

and  Weiss  (1984)  and  split  the  interval  [0,1]  into  two  parts  [0,Xg]  and  [xg,1] 
with  xg  big  enough  such  that  the  layer  terms  are  small  in  [x0,1].  We  choose 
Xg  «  2ltn  y  p(0)  «  2At0<  The  constants  in  conditions  (i)  and  (ii)  of  Theorem  3.1  can  be 
chosen  such  that  Xg  <  1  holds.  On  the  interval  [0,Xg]  we  use  T  as  independent 
variable  and  consider  instead  of  (A.1)  the  equivalent  problem 


*1  •  n  -  p  +  b1 


(A. 2) 


for  0  «  t  <  2x0, 


n,  -  *n1  +  (n(0)  +  n)^  +  b2  +  $2z  , 


P,  “  -♦P,  -  (P(0>  +  p)t,  *  b3  *  b3e  » 


•W\3 


(A. 3) 


A2*"  -  n2  -  p2  ♦  b4  , 
n2  "  **n2  +  <n  +  *ni^2  +  J  *  b5  *  b5*  » 


for  Xg  <  x  <  1  and 


Pi  -  -*,P2  -  <P  +  APl>^  -  J  +  b6  ♦  B6z  , 
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-  .** 

•  *  V  -- 


*,(0)  -  b7,  n,<0)  -  P,(0)  -  b8.  V’i  “  V  P2(1)"b10'  n2l1> 


*  b 
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(A.4)  “  VV  *  b12'  V2V  "  X*2(V  +  b13  ' 

K,(2t0)  -  n2(xQ)  +  b14,  P,(2t0)  -  P2<*0>  ♦  b,5  , 


where  *  -  tV“l'ivVV';2'3)  e  B3,X  *nd  b  "  (b1 . b15’  *  ^ 


B  -  n,»*(0,2t  )  *  (I.“<0.2T0))2  *  W1'*,(Xq#1)  *  tL"(x0,1))2  «  R  and 

3  #  ^  0  n. 


“3  \  ’  g  u 

,u  -  L-,0,2To)  x  (L’(0.2t0,)2  x  L-(*0.1>  x  <L1(x0,.))2  «  R9-  We  introduce  I  z.  3  x  on 

B  ’  end  Ibl  ,  on  B,  .  defined  as  the  au»  of  the  appropriate  norms  of  the  components 

d3,A  4,a  4, a 

Of  *  and  b  respectively.  The  bx  in  (A.2)-(A.4(  are  determined  by  the  at  intA.D 
and  obviously 

(A.S)  lbl4,l  <  U,2,X 

nolds.  An  estimate  of  the  perturbations  0j*  is  given  in 


l^mma  A.1,  Let  0S  -  (  0 . 02  *  *3*.  0 . 65* .  Sg*.  0  .  •  •  •  >  e  *4  ,  X  '  Th“" 


(A. 6) 


IBzl  .  .  <  C  Un  -r  [  +  <tn  >'b>3]«zl  x 

4*X  1  *  /b^T 


holds . 

Proof »  Comparing  (A.t)  to  <A.2),  (A. 3)  gives 

~  .  J 

02z  -  (XV1  +  X2?{  t  Xt1)n1  +  in  -  n<0)  +  Xn,  +  Xn,  )♦,  +  X  -  , 
l  ^  A  ~  .  J 

03z  -  -<XV  ♦  x2v;  +  x^)p,  -  (p  -  pCO)  ♦  Xp,  ♦  Xp,)*,  -  X  2  , 

05z  -  ({  ♦  +  xi;  «•  5^2  t  (n  +  Xn,)P>  , 

e  • 

e6 z  -  -(j  *  +  xv‘  +  v,Jp2  -  (p  +  *p,>*2  ’ 


He  use  the  Lemmas  A.d.18  and  A.d.20  to  obtain 


l’(0,2T0) 


<  t*2T0l^<l  .  +  X22T0l5;i  „  *  X2VV  -  1‘V  - 


L  <0,X0) 


L  (0,xQ) 


L  <0,»)  L  (0,2l0) 


L  (0,»)  L  ( 0 # 2T g  } 


+  ^o'Pl'  -  +  X2T  In  '  „  ]«♦,!. 

L  (0.xQ)  °  L  (0,x0)  - 

+  -j  |j|  <  const^Atn  y  / tn/b  (g—  y  +  64)  + 

+  X2tn  y  / tn/b  —  +  <S4/ In/b]  +  Xtn  y  / tn/b  ♦  6^/b  (tn/b)3^2)  + 

+  X(tn  y)2tn/b  (-  - ~  +  4^)  +  Xtn  y  /tn/b  —  +  i*/in/b)  + 


+  Xtn  4  /tn/b  ( - - —  ♦  64/b  (tn/b)5//2)  +  Xjlil.  .  < 

x  /b^“T  3' 

<  const  X  tn  t  ( - - —  +  44/b  (tn/b)3)lsl,  *  • 

*  /b“^T  3,A 


Analogously  we  get 

10. si  „  <  const  K  tn  4-  ( - - —  +  44/b  (ln/b]2)lzl,  %  . 

3  l\o,2T0)  X  /E^T  3'X 

The  Leninas  A.d.6,  A.d.19  and  A.d.20  give 


ifi.si  .  <  hj5zi  .  < 

L  (X0,D  L  <*0,1> 

<  It  .  +  *»?;■  .  +  'V  -  J'V  - 

L  (2TQ,«)  L  (0,1)  L  (2T0,«)  L  (Xq  r 1 ) 

♦  lx  -  +  -  Jx,^‘  -  « 

L  (2X0,«)  L  (2Tq,«)  L  (XQ, 1 > 

t  const[X  /tn/b  +  Xf  -1  +  54/tn/b)  +  X(^  in  y  +  S4/b  (tn/b)5//2)  * 
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+  A  -.-b  -  in  Sb  +  A[  -  — —  in  j  +  «Vb  (,ln/bj5/2 )  J  1*1 
b  "  1  /b  -  “i 


<  const  A( - - —  tn  X  +  fi4^ 

/b  -  1 


Analogously  we  obtain 

16 


2 1  <  const  A(  --■■■—  in  j  *  iVb  (in/b)5/  )***3#i  * 

L  (xrti1 )  -  1 


For  the  application  of  a  contraction  mapping  argument  to  (A.2)-(A.4)  we  analyte  the 
unperturbed  problem 


»  n  -  p  +  dj  , 


(A. 7) 


n j  -  H>n1  +  (n(0)  ♦  n)*t  +  d2  , 


?,  -  "*?,  '  <P<0)  +  *  d3  ' 


for  0  <  t  <  2T 


0' 


A2**  “  "2  -  P2  +  \ 


(A. 8) 


n'  •  ij» * n  +  (n  +  An^Jlij  +  X  +  d5  ' 


-rp  -  (p  ♦  *P,  -  2  +  d6  ' 


for  xn  <  x  <  1  and 


*,(0)  -  d7.  n,(0>  -  P,f0)  -  d0,  +2<1)  -  V  P2(,)  '  d10'  "2(1)  "  dH 


(A.9)  ’?1t2T0)  *  *2U0)  +  d12 '  *1<2V  "  **2(V  +  d13  ' 


n.<2V  »  n,(xn)  +  d 


P,<2V  “  *2l*0}  *  d  15 


where  we  introduced  the  new  inhosogeneity  d  •  (d, , . . . ,d15 )  e  ^ .  In  order  to  identify 
alow  and  fast  variables  in  (A. 8)  we  apply  the  following  transformation : 


(A. 10) 


Substitution  into  (A. 8)  yields 


*2  -  t  -  u 

"  ^2p  +  1  v 
n2  “  -(p  1)u  ♦  w 


p2  -  pu  +  w 


Au‘  *  -/ 2p  +  1  v - 


■  (2pV2£"Tl  v  +  2$‘w  ♦  3)  +  — - -  (d  -  d,)  , 

2p  +  1  2p  +  1 


(A. 11) 


Av'  -  -Ap  +  1  u  +  - ^ - -  v  +  . 1 

2<2*  +  /2p  ♦  1 


d4  ' 


Ap  +  i 


V  + 


(2p  +  1) 


2(2p  +  1) 


J  + 


Pd5  +  (p  +  Ddfi 

2p  +  1 


t* 


2*1  2J  1  ~  .  d6  "  d5 

-  v  +  - -  w  -  - J  +  -  . 

(2p  +1)  2p  +  1  2;  +  1 


As  in  Schmeiser  and  Weiss  (1984)  we  apply  a  decoupling  transformation: 


(A. 12) 


w  ■  w 


t  -  t 


*P, 


u  + 


AJ 


2p  +  1 


2/2p  +  1 


2Ap^ 

2p  +  1 


With  these  transformations  (A.7)-(A.9)  become 


%  s  .*> 


'•'6.X  on  "6,X  **  r*PUc«  'V  ,  by  X  1,1  .  The  follows, 

L  t*0‘D  1  (xo'1) 

two  UN.  give  eatinatee  for  .  in  te—  of  d  end  for  the  perturbations  E^z  in  terms 


be—  A.2i  The  estimate 


,e,6.X  4  COn,t  b^-T  ,d,4,l 


Proofi  Comparing  (A.7)-(a.9)  to  (A. 13)-( A. 16)  yields 
®i  “  ^i  for  i  ■  1,2,3  , 


'4  “25! ,  <d6  ■d5''  *5  *7=d<  ' 

/2p  +  1 


1 _  a  W  =  .  AP1  1  d5 


*6  "  - ~ -  d4  +  (p  +  - - )  —  +  (p  +  1 - — - )  - $ -  , 

2(2P  *  1)  2p  +  1  2p  +  1  2p  +  1  2p  +•  1 


2Xp  d,  -  d_  d-d 

»7  -  (1  — i!_j  e  .  > _ !n 

a5  +  1  *  1  10  *1  +  1  ' 


(pi  *  11d1Q  *  p1d11 
2p,  +  1 


_  _  j  ,  dio  "  dii 

'  12  ds  ■~2pT~rT~  ' 


e j  -  dj-l  for  i  -  8,9,13,.. .,16 


Obviously 


•el  <  const[jj|l 

L  (xQ,1) 


T= - 1  1  ,d4*  -  +  ,d5'  1  +  'd,l  .  J 

2P+1L<x0,1)  L  (x0,1)  *L1(x0,1)  6  L1(x0,1)J 


The  estimate  (A.d.  11)  implies 


4  , - - , 


2p  +  1  L  (xQ,1)  (b  +  1  )✓  1  -  x  l\ 0 , 1 ) 


b  +  1  ’ 


With  |j|  <  const(-g-^ — we  obtain 


•e,'  .  <  const  7- — — — —  I  dl 

‘^,11  b_1  4'X 


Using  the  nonnegativity  of  p  in  estimating  the  remaining  components  of  e,  the  result  of 
Lemma  A. 2  follows • 

Lewna  A « 3  ;  Let  Ez  *  ( 0 , 0 ,0  , E^z,  • . .  #E^Z/0 , 0 # 0 , z,0 , E^z r 0  # E^z,E^z)  €  B^  where 

09  « 

^  is  as  B6  \*  but  with  the  fourth  conponent  L  (Xg.1)  instead  of  L  (xQ,1).  The 


°6,X 

is  as 

B6,X' 

norm 

*e,2,x 

on  B* 

•V 

09 

r  /  —  4  V 

.  Then 

L  <x0,1) 


in  1  el 6(x  ^ 


IeIi*  ,  <  C,X( - - - -  +  iVtn/b)lzle  . 

6,A  i  Z 


(b  -  IT 

holds. 

Proof ;  Using  (A. 11)  and  Lemma  A.d.20  gives 


IE  «I  .  <  const  X(  I-  ■■  -■  -I  ^  +  lj'1  _  +  l)l  zl  . 

L  /^TTl  ‘V”  1  ‘v1’ 


<  const  X(^— — +  64/ tn/b  +  - - - r-  +  S4)lzl 


(b  -  1) 


5,X 


t  const  x(— 5— —  +  64/ tn/b)  I  zl  t  . 


IE5zI 


I*  ^ 


(b  -  1) 

<  const  X|j|l — - — - - —i  ^ 


Izl. 


j"  go  "**"  C  \ 

(2p  +  I)2  L  (XQ,1> 


<  const  X(- 


«4)iSi5(x 


(b  -  1) 


Elementary  calculations  show 


J2p 

e62  *  XU— — 


p;  v  /  2^i 

- )u  +  (- 


j2 


2Jp, 


3  -  v  .  3/0  -  5/2 

( 2p  ♦  ^ V  2p  +1  (2p  +  1 V'*  (2p  ♦  1)  ' 


)  v - - — ] —  w  +  __ — ] —  5] 

(2P  +  ir  (2P  +  ir 


and 


E  z  »  2E,z 
7  © 


From  (2.8)  we  get 


35- 


1  —  572*  1  *  "—tlS-H)  • 

(2p  +  L  (0,1) 


The  Leans  A.d.20  and  (A.d.12)  give 


,  lE?zl  1  <  const  X[(^b '  +  64/ tn/bj ( - - - -  +  64)  + 


'  V  1  *  ,E721  1 

l’(x0,1)  L  <xQ,1) 


(b  -  1)‘ 


+  — !—  +  +  S8in/b  +  ( - &1—  + 


(b  -  1) 


♦  +  i4/u^\ — s — -  +  «4)  +  ^33  ♦  «Vwe]i;ir ,  < 

b  t  V.  _  A  \2  b 


(b  -  1) 


<  const  A( - — - —  +  S4i/ In/b)  I  zl , 


(b  -  1)* 

\5  XP^ 

The  perturbations  in  the  boundary  conditions  contain  the  terms  — . . -  v  and  -  u 

2/2p  +1  2p  +  1 

at  x  -  x0  and  x  «  1.  Estimates  follow  from  Lemma  A.d.20. 

In  order  to  be  able  to  apply  a  contraction  mapping  argument  to  (A. 13)-(A. 16)  we 
analyze  the  unperturbed  problem 


(A. 17a) 


(A. 17b) 


(A. 17c) 


for  0  <  t  <  2TQ, 


*1  "  n1  -  *1  *  f1  ' 


"l  “  ,*'n1  +  (n(0)  +  n>^l  +  f2  ' 


P,  “  -^P,  "  (P<0)  +  p)*,  *  f3  , 


Au*  *  -r2 p  +  t  V  ♦  , 


Av*  *  -/2p  +  1  u  ♦  , 


V.1  .*  •*.  •*« 


v  v.' 

vv;.-/  a' 
•  si  ::vv:-:  ' 


(A. 19) 


1 


f 


( 2p  +  1 )  2 ( 2p  +  1 ) 

2J  -  1 


J  +  £6  ' 


( 2p  +  1)  2p  +  1 


J  +  f?  , 


(A. 20) 


*,(0)  -  f8,  n^O)  -  Pl(0)  -  fg,  u(1)  -  f1Q,  5(1)  -  f15,  t<1>  =  f12  , 


V2V  *  f<V  '  U<V  +  f13'  V2V  “  +  1  V(V  +  f14' 


n1(2T())  »  -(p(xQ)  +  1)u(xQ)  +  w(xQ)  +  f15  , 


*1«V  "  P<X0)u(V  +  W(X0)  +  f16  ' 


where  f  ■  (f j , . . .  ,f 16)  ia  from  Bg  ^  or  Bg  The  differential  equations  in  problem 
(A. 17)-(A.20)  contain  the  linearization  of  the  layer  equations  (2.4)  in  (A. 17),  a  standard 
linear  singularly  perturbed  system  in  (A. 18)  and  the  linearization  of  the  reduced  equations 
(2.3)  in  (A. 19).  These  systems  are  analyzed  in  the  paragraphs  a),  b)  and  c)  of  this 
appendix.  We  then  substitute  (A.a.1),  (A.b.3)  and  (A.c.1)  into  the  boundary  conditions 
(A. 20).  We  replace  tq  by  *•,  x0  by  0  and  vh(Xg)  and  uh(1)  by  0  in  the 
coefficient  matrix  of  the  resulting  linear  system  for  the  c^s.  Thus,  we  solve  a  system 
of  the  form  Ac  «  n  instead  of  (A  +  F)c  “  n,  where 

| F  L  4  c3(X/rH  '■ln^)3/2  *  exPl-^-J — )  +  Xin  X  An/b)  . 

This  estimate  follows  from  the  Lemmas  A. a. 3,  A.b.1,  Taylor's  theorem  and  the  estimates  in 
the  proof  of  Lemma  A.d.20.  The  system  with  F  »  0  is  given  by 


'  2p,  +  1  9  2(2Pl  +1)  "  ' 


i  -  +  C  —  C  -  —  t  -  r 

7  2Pl  +  1  8  9  2Pl  +  1  12  P 


=  -  t  (D  , 


(A. 21) 


1  /b  ~  1  _  „  b  -  1  b  -  1 

^  2  /  b  +  1  3  b  +  1  °4  /r*,.  . 


+  c  -  c  Kb---U  .  c 

/b(b+1,  5  7  b+1  C8 


■  fi3  -  y2v  +  2wP<xo)  -  yv ' 


ci  r "  cs  /r^r  -  f,4  -  V2V  “  yv  /Hr  ' 

b  .  _  1  /b  b  b-1 

Cl  r; -  3b+1-c4b+1°5b-1-c7b+1 

2/b2  -  1 

"  f15  "  nP(2T0)  ”  VP(X0)  b  -  1  +  wp(V  ' 


„  '  +  1  +  „  b  t  2  _  1  „  b  -  1 

1  rz -  3  b  +  1  C4  /r-,.  ...  5b-  1  C7  b  +  1 


2/b  -  1 


/b(b  +  1) 


f,,  “  P  ( 2t  )  +  V  (x  )  — - r  +  w  (x.) 

16  p  0  p  0  b  -  1  pO 


Lengthy  calculations  show  that  the  coefficient  matrix  satisfies 

|a'1l<c4F4t. 

With  the  constants  in  Theorem  3.1  chosen  aPProPriately,  the  estimate 

holds,  which  implies  the  nonsingularity  of  A  +  F  and 

l(A+F)"tL<2C4rAT 

Using  the  estimates  on  the  particular  solutions  in  lemma  A.a.1,  in  (A,b.4)  and  (A.c.2)  we 
obtain  an  estimate  of  the  righthand  side  n  of  (A. 21): 


‘v'W’V 


. .N 


-V '.'•V-  .  - 


I'll.  <  con»t['b'  j  -  +  tn/bjlflg  ^ 

This  implies  for  c  “  (ci,...,c9): 

(A.22)  |c|  <  const ( - - - 5  ♦  tn/blifl.  . 

(b  -  I)2  6'X 

Using  (A.22),  tbs  estimates  (A.b.4),  (A.c.2)  and  Leans  A. a.)  we  obtain 


ill.  .  <  const( - - - r  +  (tn/b)2)lfl,  . 

5,X  (b  -  I)3 


The  difference  between  1*1^  ^  and  I • I  *  ^  only  appears  in  the  fourth  component.  Using 
(A.b.S)  we  immediately  obtain 


(A. 23) 


1*1. 


5  X  4  C5l - 3 - 3  +  \ 

S  (b  -  I)3  6,A 

A. 3  and  (A. 23)  imply  the  condition 


(A. 24) 


C  C  X( - - -  +  S4(ln/b)5/2)  <  ^ 

(b  -  1) 


for  the  applicability  of  a  contraction  mapping  argument,  with  the  constants  in  Theorem  3.1 
chosen  appropriately  (A. 24)  is  satisfied  and  the  solvability  of  problem  (A. 13)-(A. 16) 
follows.  Furthermore  we  obtain  the  estimate 


and  from  Lemma  A. 2 
(A. 25) 


1*1.  ,  <  const( - - - -  +  (»n/b)2)lel,  .  , 

5,X  (b  -  I)3  6'X 


1*1.  ,  <  const! - - - 7  +  (in/F)2)ldl  .  . 

5'A  <b-1>4  4'X 


The  transformations  (A. 10)  and  (A. 12)  give 


Thus,  (A. 25)  implies 
(A. 26) 

If  the  condition 


,Z,3,X  4  COn8t  iT7  ,Z,5,X  • 


«*•-  l  4  C.( - ' - k  +  <*n/b)2)ldl  . 

3<X  b  (b  -  1 ) 3  4'X 


ciVln  f  ^-."11/2 +  *<*"**>*) 4  i 
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is  satisfied  (which  again  is  guaranteed  by  the  assumptions  in  Theorem  3.1)  we  can  apply  a 
contraction  mapping  argument  to  obtain  solvability  of  problem  (A.2)-(A.4)  and  the  estimate 


III,  .  <  const( - - - jr  +  (tn/b)2)lbl.  ,  . 

3'  (b  -  I)5  '* 


Thus,  problem  (A.1)  also  has  a  unique  solution  y.  The  differential  equations  in  (A.1)  can 
be  used  to  show 


and  the  estimate 


(A. 28) 


"y*  1.A  4  Const  b^T  ,Z,3,X  ‘ 


Finally  (A. 5),  (A. 27)  and  (A. 28)  imply  the  validity  of  the 
Lemma  A. 4 :  Problem  A.1  has  a  unique  solution  y  and  the  estimate 

lyl,  .  <  const( - 1 - r  +  <tn/b)2)lal,  . 

1,A  <b  -  ir  2,A 


a)  The  linearised  layer  equations. 

We  consider  (A. 17).  Integration  of  (A.17b,c)  and  using  (2.6),  (2.8)  yield 


~  e  b  ~ 

"l  “  C3fi  +  8  *1  +  "q  ' 


(A. a.1) 


c4  1  -i>~  .  ~ 

P1"^“  "  —  *  *1+P*' 


where  C3,  are  constants  and 


n  (X,  -  /  e+(T)-*<s)f,(s)ds  , 

q  2To 


A — 


Pq(T)  -  J  e-*(t,+*<B,f3(s)ds 


The  monotonicity  of  (1  immediately  implies  the  estimates 


(A. a. 2) 


'V  -  «  **2*  1  '  ‘V  -  «  »V  1 

4  L  (0,2tq)  L  (0,2X0)  q  L  (0,2T0)  L  1  ( 0 , 2Xg ) 


Substituting  (A. a. 1)  into  (A.17a>  givas 


(A. a. 3) 


•T  ♦ 


4_  .  v*  -  *  S,  -  i, . ,, 


The  homogeneous  part  of  (A. a. 3)  has  the  linearly  independent  solutions  end 

*  T  *-2 

F(t)  J  ♦  (s)ds.  For  the  construction  of  a  particular  solution  of  the  inhomogeneous 
0 

equation  we  use  the  results  of  Fife  (1974)  and  get 


t  „ 


-  c^{2x0)1Hr)  j  <f2(s)ds  ♦  c2F(D^_1(0) 


-  *(T)  J  *"2(s){J  ♦(u)[c  e*(u)  -  —  e_,il(u)Jdu 

0  a  /b 


2X. 


♦  /  4»(u)[n  (u)  -  p  (u)  ♦  f1(u)]du}ds 


(A.a.4a) 

-  c  *(2t  >»(T)  ;  *"2(s)ds  ♦  c  *<t)*_1<0> 

0 

-  c  i(t )  /  (1  -  e*(,))*-2( s)ds 

J  o 

’  T 

A  -*(«)  1  *-2  ~ 

-  c.*(T)  J  (e  '  -  1)  —  *  (s)ds  +  *  (T)  , 

4  0  /b  p 

~  *  T  "_2  .2*0  > 

where  ♦  (T)  **  J  ip  (s)  j  ♦<u)  [n  (u)  -  p  (u)  +  f  (u)]duds.  Differentiation  of 

P  0  s  q  q 

(A. a. 4a)  gives 


(A. a. 4b)  ♦)(T)  -  Cj<i(2Tq  )  [  (n(T  )  -  p(T))  j  *  2(s)ds  +  ^  ^T)]  +  c2(r(T)  -  p(T  )  )4>- 1  (0 ) 


*-1 . 


A- 1 
\,U  »  / 


cJ(n(T)  -  pit))  j  (1  -  el|'(s>)t"2(8)ds  +  (1  -  e*11 5  J*'1  (r  )  j  - 
0 


41- 


-  c4l<n<T)  -  p(T)>  j  (e-*  -  1  >  — 1  *‘2(s)da  ♦  <e‘*(T> 

0  /b 


1)  —  *“1{T) 
✓b 


+  VT>  . 


The  estimates  of  the  solution  are  collected  in 
Lemma  A. a.  1: 


a)  'V  -  <  const(/wi  |c  |  +  Jc  |  ♦  Wb  |c,|  ♦  Wb  |c.|  ♦  in/b  lflc  ,) 

L  (0,2t^)  J  4  6, A 

b>  I*/  „  <  const( ( 1  +  /in/b)  |c  |  +  |c,|  +  in/b  |c.| 

L  <0,2t0)  1  i  b  -  1  2'  '3' 

+  in/b  |c4|  +  (1  +  in/b)lf lg  x)  , 


\c2\  + 


c)  In  I  m  <  const ( ( - - —  +  /in/b)  |c,|  +  — b  4  |< 

L  (0,2tQ )  /b"^T  1  b  ‘  1 

♦  (1  +  in/b)  jc3|  +  (1  +  in/b)  |c4|  +  (7  +  in/b)  Iflg  A)  , 

d)  ,?1,L-(0,2T  )  '  con8t(/r?r  lcil  +  b^T  ic2l +  hi  +  lcJ  *  • 


Proof :  The  results  are  obtained  by  using  the  Lemmas  A.d.4, 
and  A.d.12  and  the  estimates  (A.a.2). 

With  np(T)  -  e*(T,*p(T)  +  nq(T)  and  pp<T)  -  - 

define 


Aede5#  AedeBp  Aed>  10|  A  «  <1  e  1  1 


e-*(TV>  *p  ,T) 

p  q 


W 


VT)  -  ?,(t)  -  yT>  , 

r^tx)  -  nj  (T)  -  np(T)  , 
Ph<T>  “  P,(T>  -  pp(T)  . 


We  have 


D 


li®  \(T)  ■'u/M!-c3b+,tC</S(b+1) 


b)  li®  *h<t)  -  c,  -j  , 


c)  li®  ^<T>  -  c,  -—rr  +  cj  g-rr  -  c4  b^rr  ' 


t+« 


2/?T7 


d)  liB~h(t)..Ci_^+C3^  +  C4_^_ 


,’s*  / 

/*  .  r  e  "JL 
1  •  •  4  •  a  ' % 


VT;  ■  J 


Proof:  The  Leonas  A.d.3,  A.d.5,  A.d.8,  A.d.9,  A.d.11  and  A.d.12  innedlately  imply  a)-d) 
A. a. 3: 


a)  l*h(2t0)  "  *h(">l  *  con8t  X  /g T~T  (tra^l3/2 


b>  l?h<2T0)  "  ?ht“>l  *  con8t 


c)  l“h<2T0>  "  "h<“)l  <  °°nst  X  /g"V~|  Un/b>3</2 


dl  IPh<2To)  "  Ph(-,l  *  conBt  x  1  Un/b)3/2 


/7T7 


Proof :  The  results  can  be  shown  by  combining  the  Lemnas  A.d.6,  A.d.9  and  A.d.11- 17. 
b)  Analysis  of  (A. 18) 

The  analysis  of  (A. 18)  is  facilitated  by  the  transf ornation 
(A.b.1)  u  -  u  t  v  , 

v  -  U  -  v  . 

Using  (A.b.1)  we  obtain  fro®  (A. 18): 


Ap •  “  -/ 2p  +  1  U  +  j  (f^  +  tg), 

(A.b.2)  _ 

IV  -  / 2p  +  1  V  +  y  (f  -  f.)  . 


The  behavior  of  solutions  of  (A.b.2)  is  characterized  by 


Le— m  A.b.1;  Let  ly'  •  -o(x)y  +  f(x),  x  e  [a,b],  a(x)  >  a  >  0.  Then  there  Is  a 
solution  y^  of  the  homogeneous  equation  and  a  particular  solution  y^  which  satisfy 

a 

■y,,1  „  <  '»  yh(a)  “  1#  lyh{b,l  *  exp("  x  (b  ”  a))  * 

L  (a,b) 

'V  .  4  i  •*»  -  *  *v  -  «  *_1*«  , 

p  L  (a,b)  -  L  (a,b)  p  L  (a,b)  L  (a,b) 

-1  *  x  x 

Proof:  We  use  ( x )  «  exp(-A  j  a(t)dt)  and  y  (x)  “1  /  exp(-A  j  a(s)dsjf(t)dt 

a  p  a  t 

An  analagous  result  holds  for  a(x)  <  -a  <  0.  Thus,  we  have 


(A.b.3 ) 

u  -  c5)ih  +  c6vh  +  pp  +  vp  . 

v  *  cSwh  -  Vh  +  vp  ~  vp  * 

where 

(A.b.4) 

lul  ,  Ivl  <  const (|c_  |  +  | c,  |  +  Ifl,  , 

L  (x0,1)  L  (x0,1)  5  b  6'A 

and 

(A.b.5) 

lul  „  ,  <vl  „  <  constt  |  c,  |  +  |c,  |  +  Ifl*  .)  . 

L  (xQ,1)  L  <x0,1)  S  6  *'* 

c)  The  linearized  reduced  equations 
Integration  of  (A. 19)  gives 


=7  - - °g 

2p  +  1 


+  w 


2(2p  +  1) 


(A.c.1) 


t  »  c, 


7  2p  +  1 


+  C8  '  c9 


+  t 


2p  +  1 


where 


“  1  *  • 

w  (x)  - -  /  (2p(t )  +  1)f  (t)dt  , 

P  2p(  x)  +11  6 


t  (x)  -  2w  (x)  *  }  -  2t, (t) )dt 

0 

Since  (2p  +  1)  la  aonotonically  decreasing  the  estimates 

Iwl  m  <  conat(|c7|  ♦  |c9|  +  Ifl^  x)  , 

X.  (x0«D  ' 

(A.c.2)  Itl  m  <  conat(|c7|  +  |c8j  +  jc9(  +  I f I 6  x)  , 

L  (xQ,1) 

|J|  -  |c,| 

hold. 

d)  Collection  of  technical  results 

He  now  state  sone  useful  general  results* 

* »  A.d.1:  Let  y  be  the  solution  of  the  initial  value  problem 

y  -  a(T)y  +  f(x),  x  >  0,  y(0)  -  yQ  , 

with  S/f  e  CIO,-)*1,  a(x)  <  -u  <  0  for  x  >  0.  Let  y.  -  lim  ~y.  Then 

y+- 

|y(x)  -  yj  <  |y  -  y-|e”WT  +  i  e“MT/2lf  +  ayj  „  ♦ 

(A.d.1)  L  ‘°'T/2) 

+  i  ,f  +  •*.*  • 

*  L  (X/2,X> 


holds. 

Proof  >  He  set  u  »  y  -  y„.  Then  u  is  the  solution  of 

u  -  a(x)u  +  h(x),  u(0)  -  yQ  -  y,,  , 

where  h(x)  <*  f(x)  +  y^alx).  Thus, 

X  XX 

u  *  (yfl  “  y„)exp(J  a(s)ds)  +  /  exp (/  a(u)du)h( s)ds 
0  Os 

He  split  the  integral  in  the  particular  solution  into  two  parts  and  use  the  upper  bound  of 
a(x)  to  obtain  the  estimate 


)  g  e  C[0,“J  <»■>  (g  6  C[0,«)  and  lim  g(x)  exists  and  is  finite.) 


which  implies  (A.d.1). 

Lemma  A.d.2i  Let  y  be  the  solution  of 


y  «  a<T)y  +  f(T),  T  >  0,  y(»)  -  y  , 


a,f  e  Ct0,»],  a(x )  »  u  >  0  for  t  >  0  and  y  "  -lim  - ,3-r. 

t*~  *(T' 


(A.d.2) 


|y(t)  -  y„l  <  £  ♦  «yj 


L  (T,«) 


Proof :  The  function  u  «  y  -  ym  is  the  solution  of 

u  -  a(T)u  +  h{T),  u<“)  «  0  , 


|ult)|  <  |yQ  -  y„|e_WT  +  J  e“u(T_8) |h(s) |ds  < 

0 

«  |y0  -  y.le~l,T  -•»>*«  i  ‘••UT/2  -  •-ut>  ♦ 

L  (0,t/2 )  " 


L  (T/2,1) 


I  (1  -  .“UT/2  )  # 

P 


where  h(T)  »  f(T)  +  yaa(T).  Thus 


t  T  *  s 

u  -  J  exp(/  a(u)du)h(s)ds  -  -  j  exp(~  J  a(u)du)h(s)ds 
•  S  T  T 


Using  the  lower  bound  of  a(T)  we  obtain  the  estimate 


|  u(  t  )  |  <  j  e  U<8  T)|h(s)|ds  <  •j-  Ihl  m 


L  (T,») 


Lemma  A.d.3;  Let  g  e  C[0,«);  f,  h,  f(T)  /  g(s)ds  e  C[0,»J,  and  f(-)  -  0. 

0 

Then 


%v. 


.**  /•  .**  a*V**  •* 


(A.d.3) 


lim  f(T)  J  g( 8 )h( b )da  -  Urn  h(T)f(T)  J  g(s)ds 
T+“  0  T*»  0 


holds. 


Proof  t  The  limit  on  the  right  hand  aide  of  (A.d.3)  exists.  It  remains  to  show 


T  T 

|y(T)|  -  |h(t)f(T)  J  g(s)da  -  f(T)  /  g(s)h(s)da|  <  e  for  t  >  t(c) 
0  0 


He  set  Ihl  „  «  M.,  If ( T )  J  g(s)dsl  K  “  M-.  Clearly 

L  <0.«)  0  L  (0,-) 


|y(t)| 


J 

0 


g(s) (h(T  )  -  h(a)  )ds | 


t 


t  T 

<  |f(T)  J  g(s)(h(T)  -  h(s) )ds|  +  |f(T)  £  g(s)(h(x)  -  h(a))| 
0  t 


holds.  We  choose  t  “  t(e)  so  that 


|h(t)  -  h(s)|  <  for  t,s  >  t 


holds.  Then  we  choose  t  »  t(e)  >  t  such  that 


|f(x)|  <  e( 4H.  J  Ig(s)jds)"1  for  t  > 
0 


This  implies 


| y ( t  J |  <  e  for  t  >  t  . 

The  following  Lemmas  are  results  on  the  behavior  of  the  layer  solutions. 
Lemma  A.d.4:  There  are  positive  constants  Cj  <  c2  such  that 


,/in/b  * 


|>(0)  «  c, 


J tn/b 


Proof :  The  relations  (2.11)  imply 


i(o> 


♦  in/b  - 


b  +  1^ 
b  -  1J 


1 2(tn/F  -  ^ - -)  — 

/b  +  1 


/ 


2  in/b  F(/b), 


with  F(5 ) 


e  -  i 

(€  *  UinC 


Elementary  calculations  give 

f  e  c[i,-l. 

This  implies  the  existence  of  Ci,c2 


F(£)  >  c  >  0  for 
>  0  such  that 


2 

C. 

—  <  F(/b)  < 


5  e  ii,-) 


Lemna  A.d.5; 


_  <  J(o), 

n  -  p  L  ((>,«•) 


T-*»  n  -  p 


/ 


b  -  1 
b  +  1 


Proofs  Let  H(^) 


-a)2.  With  (2.11)  we  get 
n  -  p 


H(-in/b)  -  *(0)2 


To  verify  the  second  assertion  of  the  lemma  we  use  the  rule  of  de  l'Hopital  and  the 
relations 


( A.d.4) 


di  n  -  p  dn 
-  -  — — sr 

«♦  J  d* 


+  n. 


£ 

di> 


1 


which  follow  from  (2.11)  s 


* 

lira  H(<|i) 
♦♦0 


lira  — 

<|/*0  (n 


lira 

♦♦0 


_ 2(n  -  p) _ 

2<n  -;>(|4r  +  "  +P) 


b  -  1 
b  +  1 


To  complete  the  proof  of  Lemma  A.d.5  we  show  that  H< ♦ )  is  monotonically  decreasing: 


hold  for  1  >  0. 

Proof ;  z  «  ♦  la  the  solution  of 

a  a  • 

A 

z  -  —  ~-£  z,  z(0)  -  <.(0) 

* 

*  -1 

An  application  of  Lemma  A.d.1  with  V  *  H»(0)  (see  Lemma  A»d.5)  yields 

(A.d.5)  ♦<!)  <  t|/(0)exp<  -T/i(i(0) )  . 

Taylor's  theorem.  Lemma  A,d.5  and  (A.d.4)  imply 

|+<T)|  <  *<0)*(T) 

|n(T)|  <  const  ♦  (<>  )♦(*  ) 

|p(T)|  <  conBt  g— — p  l|l(0)lt>(T  ) 


This  completes  the  proof  of  Lemma  A.d.6 


Lgga  A.d.7: 

1> 


<  const  X  <>(0)  holds  for  1  >  0. 


L  <iV-) 


Proof :  We  set  “  (—  ~  ?) 


n  -  p,2  (F'<*)) 


-«  where  P(  P> 


An  application  of  Taylor’s  theorem  yields 


2F<*) 


b  *  .  1  -<>  :  b  ♦  1 

— — r  «  ♦  - - r  e  T  -  *  - 


b  -  1  b  -  1 


which  implies 


Hence 


»<♦>  *i2r^ri1  +  • 

H<i)  -|44[,  +  otiH4)J  • 


The  assertion  follows  from  Lemma  A.d.6. 

:  t  ; 

Lemma  A.d.8;  Let  y  -  <Kt)  j  ♦(s)  ds.  Then 
~  0 

•>  u"  *(T)  -  /i¥l 

X+m 


b)  lyl  .  <  *(0) 

L  (0,-) 


Proof :  a)  de  l'Hopital's  theorem  and  Lemma  A.d.S  give: 


lim  y(T)  *  lim  (-  n— ~  P-J  — 

T~  J 


-  lim 

T+®  n  -  p 


_  »  -  /b  ~  1 
11  /  b  +  1 


b)  y  is  the  solution  of 


y  -  n—  p  y  +  i,  y(0)  -  o 


b  -  1 


Since  y  “  0  holds  at  a  relative  minimum  or  maximum  of  y  in  (0,®)  we  obtain 


U— wi  <  ^(o ) 

n  -  p 


From  y(0)  =•  0  and  a)  we  conclude  b). 


*2  1  A» ^ 

Lemma  A.d.9:  Let  y  «  ♦  (t )  /  t  (side.  Then 

0 


a)  Urn  y(t)  -  I  /fZI 


b)  iy«  „  <  *(0) 

L  <0,«) 


c)  «y  -  1  /^-il 


2  ✓  b  +  1  ® 

L  (2it0,«) 


<  const  1 0 >  for  i  >  0. 


Proof:  a)  From  de  l'Hopltal's  theorem  and  Lemma  A.d.5  we  derive: 


Urn  y(T>  -  Urn  (-2 

X-**  ,3  ,2 

*  V 


b)  +  is  positive  and  monotonically  decreasing.  Thus,,  Lemma  A.d. 


»  T  *  -  - 

y(T)  <  ♦  (T  )  /  $  ( s ) ds  <  \K0)  . 

0 


c)  y  is  the  solution  of 


y  “  2  -n-~  ^  y  +  1,  y(  0 )  *  0  . 


An  application  of  Lemma  A.d.1  with  W  =  — -  (see  Lemma  A.d.5)  c  ves 

*(0) 


8  gives 


-51- 


<fr(0) 


L  (0,-) 


*<0) 


*  A  _ _ _ 

1  *(0)11  +  A._-J| 

2  x  /  b  +  1 


L  (T/2,-) 


which  implies 


ly  -  1 

Y  2  /  b  +  1 


<  const 


L  (2it0<-) 


X<1  +  itO)X2i  +  *(0)  /p-~  i2(0)XiJ 


A.d.10^.  Let  y(T)  -  -J —  /  *(s)ds.  Then 


+<T) 


■y1  *»  <  const  /l n/b 

L  (0,-) 


Proof:  Obviously  y(x,  -  -*(t,/*(t)  and  y  le  the  solution  of 


y  "  ■  y  -  i.  y(0)  -  -t<o)/i<o) 


y  assumes  its  maximum  at  t  -  0  or  at  a  stationary  point.  Using  Lemma  A.d.4  w.  obtain 


y(0)  ■  <  const  / tn/b 


*(0) 


For  stationary  points 


|y(T)|  -  l^-^l  <  *(0 )  <  const  4n/F 
n  -  p 


holds  by  Lemma  A.d.5 


A.d,  1 1 1  Let  y  *  — ~ 


a)  lim  y(T)  -  /±I1, 


b)  lyl  —  <  const  /ltv/b, 

L  (0,-) 


c)  ly  -  /g  -~-|l  „  <  const  iVlO)  for  i  >  0. 


L  (iT0<-) 


Proof i  a)  We  apply  de  l'Hopital's  theorem  and  Lemma  A.d.5s 


lim  y(T)  -  11m 


-  -  lim  er 

^♦0  n  -  p 


b)  y  is  the  solution  of 


y  -  -  y  -  «*’.  y(°*>  “  /I'H' 


As  in  the  proof  of  Lemma  A.d.10  we  have  to  estimate  |y(T )  j  at  t  »  0  and  at  stationary 


points ; 


r~ 

y<0)  -  (l  -  1(0)  <  const  — - V. 

'/b  /b/in/b 


Similarly  to  the  proof  of  Lemma  A.d.4  we  obtain 


y(0)  <  const  /in/b  . 


For  stationary  points  we  have 


J y |  »  e1*1 1  >'  —  irj  <  V(0)  <  const  /£n/b  . 


c )  We  apply  Lemma  A.d.2  with  P  “  P  fO)  and  obtain 


ly(T>  -/rrrl  <  “ 6,1,1 


L  (T,») 


which  implies 


ly  - 


fb  - 


4  *-7* +  /i^H1  -  %  +  *<0>M  -  •*'  - 


b  +  1  L-UV-,  '  ■  “  •  L  (it0,«)  L  (iT0,«) 


<  const  ^i(O) 

by  the  Lemmas  A.d.6  and  A.d.7. 


[/%T~J  ^>2(0)  +  Xi*2C0)  J 


Lemma  A.d.12:  Let  y 


-  1 


__  ft 

✓b  ♦ 


Then 


a)  lim  y(T)  - 


J.-.P-.Z. 

/  b(b  + 


-  1 


b)  lyl  . 

L  (0,«) 


b(b  +  1) 

<  const  /tn/b 


c)  ,y  “  / b(b  +  1) 


,113 


<  const  1  i|i  (0)  for  i  >  0. 


-ill  ill 

e  T  i  -  e 

Proof;  a)  y(T)  «  -  •  — ; - .  Lemma  A.d.lla)  yields  the  result. 

fb  1 


(I) 

b)  follows  from  Lemma  A.d.llb)  and  from  e  <  /b. 


<  *  +  l»(0)A  const  /in/b  + 


+  itoji  -  azl  ♦  *«»  i /im  + 

*  L  (1V  >  J  t-  ^v-) 


const  Xi[Xi  g  ~  ]  +  *(0)  /tn/b  +  *4<0)  +  |  ~  ■■.1  J3(0)]  . 


*  T  -*(»)  _ 

Lemma  A.d.14i  Let  y(T)  *  x )  J  - - 

0  /b^s, 


ds •  Then 


a)  lia  y(T)  -  _£ - ! — , 

T*»  /b(b  t) 


b)  ly  - 


/b(b  +  1)  L* (21t0,-) 


<  const  XH4(0)  for  i  >  0. 


Tb®  proof  of  this  Leaaa  proceeds  analogously  to  that  of  Leaaa  A.d.  13  and  is  therefore 
oaitted. 

*  ^  ^  ^  * 

L®—  A.d.  15 :  bet  y<T)  -*(T)(n(t)  -  p(T))  /  ♦"2ls)ds.  Then 

,  0 
a)  lia  y(T)  -  -  -I 

T+m  * 


b)  ly  ♦  -j*  „  <  const  lii3(0)  for  i  >  0. 


L  (21T0,-) 


n  -  p  *2  T  *-2 

Proof:  a)  y  «  ♦  /  *  (s)ds.  He  coabine  the  Leaaas  A.d.5  and  A.d. 9a)  to 

i  0 

obtain  a ) .  T 


1  ,  .  n  -  p  1  /b  "  1  ■  ♦ 

b)  *y  +  j1  •  <  ly  -  2  /  b  +  1  -,,4T 

2  L  <21Tn.-)  ;  L  <2iV  ’ 


.  ,n  -  p  1.  /b  ~  Jl  +  1,  < 

1  2/b+1  2  L*<2iTA,-> 


*  -  ;  t  •  ,  /jtt 

<  |^-£|  .  »♦  «*>)  *  <s)d*  ■  2  / FT11 

*  L  (2iT  ,«)  0  L  t2iTo'  ' 

^  0 


*  1  /b  -  1|  n  -  p  +  /b  *  1| 
l/b+1  *  /b-1  j1-(2it  ,•) 


<  «...  •  /f-H 


' ,  -  .«•> 


L.dt.te.  Let  y(T>  -  (|.<T)  -  pit))  /  — — - dB*  Then 


0  *2<s> 


/b-1 

a)  U»  y(T)  -  -/ b  — 

T-mp 


b)  »y  «■ 


/b  -?i|  <  const  Xi/|^4  +4(0)  for  1  >  °* 

/___•  /b-1 


*  *  1  e"*l“}  -  1 

k.d.17S  Let  y(t)  -  (n<t)  -  p(T))  /  - j -  ds.  Then 


0  /£•*{■) 


a)  U»  y(T)  -  -/b(£"+  IT' 

,y  /b(b+1)  L“(2iT  ,•) 


const  **<°>  for  1  >  °* 


Th«  proof >  of  tht  L®RM8  Aede16#  A*d*l7  srfl 


similar  to  that  of  Lettw&  A-d-15 


A-d-li8  *>  <  const  (^-U-  ♦  64/b(i„^,5/2) 


b)  'V  .  <  conat(^  ♦  64/b(tnVb)5/2), 

L  (0,-)  b  J 


c)  In^  w 


<  con«t( - —  +  «Vb(in/b)5/2), 

t  (o,«)  /T^T 


d)  ,pl'  ■  <  con»t( - - —  +  a4/b(tn/b)3/2h 

t  (o,-)  /bTT 


Proofs  The  estimate. 

(A.d.6) 

A  A  *■  ■ 

,n*  i  •  •P1  .  <  const/Wb 

!•  (0,-)  1  (0,») 

(A.d.7) 

1  <  const (in/b)3//2. 

L  (0,-) 


hold,  since 


Ini 


1  *  ,pl  i  <  In  -  pi  •  *(o) 

L  <0,«)  l\o,-)  L’tO..)  ^ 

hold.,  (A.d.6)  follow..  Using  (A.d.5)  in  the  proof  of  Le».  *.d.6  w.  get 


,T*1  ,  ^  <  *(0)IT  exp(-T/*(0))l  1  .  J3(0) 


L  (0'“*  L'(0,») 

which  implies  (A.d.7).  Now  we  are  able  to  find  estimate* 


on  the  functions  In  (2.16) 


sing  P  (0)  ^(b  _  1}2  "  5  ^  b  ♦  1  and  ***  monotonicity  and  sign  of  we  get 


In  I 


*  const  V—  +  64<£n/b)3/2} 
/.  * 


,  <  const  ( - - - j  ♦  «*)  K4  *  +  <*«^>3/2)  « 

P  L* ( 0 »*• )  (b  -  D2 


<  con8t  ( - ?—  ♦  «4/b  Un/b)3/2) 

✓b  -  1 


„e  set  H<t>  -^COM.*"1  -  e-*<t))  +  np(x>  -pp<t).  use 


•  (0)  ..(__b_^-S4)(|^)^(0) 

1  (b  -  1) 


and  obtain 

IHI 


.  a-  ,  b  .  s4wb_=_li2yfHl/C  *  - 3 —  +  «VbUn/b)3/2| 


<  const  (  ■  J—  +  «4/b  (in/b)3/2)  • 
/b  -  1 


This  estimate  and  the  Lemmas  A.d.8  and  A.d-10  imply 


•*  * 


<  const  ln/b  f - *  54/b  <in/b>3/2)  < 


p  l“<o,«)  ^b  -  1 

v  const  ♦  «Vb  Un/b)5/2) 


Obviously 


1  .  .~JP.  J)  -  *_1(T)  /  +(s)H(S)ds  , 

VP  *  P  X 

* 


which  implies 


.  .  <  const[/ffI  ♦  64/S  (tn/b,5/2)  * 

P  L  <0,-) 


♦  /!^  (  1—  +  i4/b  <in/b)3/2)]  < 

✓b  -  1 


<  constfr-  +  i4/b  (tn/b)  ^  ) 


The  constant  c  in  (2*15)  satisfies  (see  2.1) 


=  1  ■  |2( - ~ - J  -  <4Hg-+  |)2»<0)|  <  oon»t(^~— ♦  64/tn/b) 

(b  -  1)  ° 


Thus,  we  have 


'V  -  «  const(I*LI— 1  ♦  fi4/tn/b  +  £EZj.  +  ««/£  ,tn/b)5^2)  < 


L  <0,«) 


<  conat(- b  ~  1  +  fiVb  <Wb)5/2) 


'♦t',*.  „  “  ,CV-E+  V  ...  .  <  const((£ZI  ♦  S4^  _L_  + 


L  (0,-) 


*<0) 


L  (0,-) 


+  j"  +  4 tn/b^ ) 5/2 )  <  const  (-^  +  64/b  Un/b)5/2J  , 

,ni'L-(o  -j  <  «4/b(in/b)S/2) 

*  (-i-j  *  ♦  _J_  .  i‘,,^,3/2]  < 

(b  -  1)  /b  -  1 

<  const( - 1—-  +  «4/b  (tn/b)5/2)  , 

/b  -  1 

,Pi'l“(0  -)  *  con"t[b^_T  (^T"1  +  fi4/b  <*n/b>5/2) 

+  /b  (— b  ♦  «4)(|^4)2/^  ♦  -L-  +  54/b  (tn/b)2/2] 

(b  -  1)  /b  -  1 

<  const ( - - —  +  «Vb  (in/b)3/2l  . 


jMgm*  A.d. 19i  •)  •  ♦.I  m  <  const  Xi(in  ■  +  S4/b  (tn/b)5^2) 

L  (2it0,-) 

b)  »♦,!  <  const  X1^  tn  y  «■  «Vb(ln/b)5/2)  , 

L  <2iT0,-) 

c)  In.l  <  const  X1( - 1 —  tn  y  ♦  «Vb(tn/b)5/2)  , 

1  L  (2iT0,«)  •'b  -  1 

d)  Ip.!  <  const  X1  -  ■  -  tn  y  . 

1  L  (2iT0,«)  /b  -  1 


hold  for  i  >  0. 

Proof i  The  estimates 


(A.d.8)  Ini  ,  Ipl  ,  <  const  xVtn/b 

L  (1T0,«)  L 

(A.d.9)  It*I  <  const  X^n  y  (tn/b)3//2 

I.  (1T0,«> 


hold.  The  proofs  of  (A.d.8)  and  (A.d.9)  are  similar  to  those  of  (A.d.6>  and  (A.d.7 
therefore  oaltted.  Proceeding  along  the  lines  of  the  proof  of  Lemma  A.d.18  gives 


In  I  . 

P  L  (IT,,,-) 


V  - 

P  l  (1T0.-) 


<  const  xHn  -f  ( - * —  +  fi4(tn/b>3/2)  , 

*  /b^T 

<  const  X^tn  y  +  64(t n/b)3/2)  , 


IHI  m  <  const  X^tn  y  (* —  +  64(  tn/b)3//2j  . 

L  UTn,«)  /b  -  1 


Mow  we  apply  Leona  A.d.10  and  obtain 


•v 


<  const  /ln/b  lyl 


)  and 


(A.d.10) 


L  (2iTn,«) 


L  (2iT„,-) 


ds .  y  is  the  solution  of 


where  y(t)  •  *(x)  J  *‘1(s)lHl  „ 

0  L  <s,») 


y  .  £  y  +  IHI  m  ,  y(  0)  «  0  . 

I  L  <s,*> 


*•1 

An  application  of  Lemma  A.d.1  with  M  *  ♦  (0)  gives 


|y(T)|  <  i>(0)  exp(-i_1(0)  yJiHi  „  +  _ 

L  (0,-)  L  (1/2,*) 


Using  this  estimate  in  (A.d.10)  implies 


I*  I 

n  «* 

P  L  (2iTQ ,“) 


I*  '  . 

P  L  (2iTQ,-) 


<  const  X^in/b  [ - —  +  64/b(^n/b)3/,2  +  tn  -J-  ( - 2 —  +  S^in/b)3^2)  | 

/b  -  1  /b  -  1 


<  const  Xi[ln  -y  ^  -*  +  64/b  (ln/b)5/,2j 


<  const  X1,  ^ *n  X  b  1  +  64/b  (tn/b)5^2)  + 


+  J In/b  Xitn  Y  (  -- -  +  54(in/b)3^2 )  ]  < 

/b  -  1 


<  const  X^i  *n  X  +  <tn/b)5^2) 


These  estimates  enable  us  to  show  the  assertions  of  Lemma  A.d.19: 

It  I  <  const [ *  *  64/ In/bJ  — * —  X2^i(i(0)  + 

L“(2iV°,) 

♦  X^tn  I  ^3  +  S4/b  <ln/b>5/2)]  , 

I*  I  „  <  const [ (-r  +  «4)X2i  t— — — r  J2(0)  +  X4(I  *n  y  +  «4/b  Un/b)5/2)]  , 

L  (2iTQ,-) 


e4t  >‘1*”  ^  ' 

L  (2it0.-) 

L  <2iT0,») 

the  lut  Lemma  of  this  paragraph  gives  estimates  on  the  smooth  terms  in  the  asymptotic 
expansion. 

Lemma  A.d.20;  Let  0  <  x  <  1.  Then 


•>  '5".  b’  ,iV, *  ““‘l 


L  (0,x) 


c)  lp"«  .  .  <  con»t(irf-T  +  <  )» 

L  (0,x) 


d)  »P*«  .  * 

L  (0,1) 


<  const(-^-g—  ♦  «Vt„/b),  f)  i} 


d)  «pi»  .  . 

L  (0 ,x) 


g>  ^'l-<o  i)  *  C°n8^ — ^~+  6  An/*  *  h> 


2 

(b  -  1) 

+  0 

[  b2 

+  t 

l  4 

(b  -  1) 

r  1 ...  * 

6\ 

✓b  -  i 

(  1  m 

To 

(b  -  1) 


Proof;  We  will  use  the  estimates 


( A.d. 11) 


1“ — - I  « - 

2p(x)  +1  (b  +  1)/t  -  * 


i  .  _ _ b  -  1  | _ 1 _ I  <  const  , 

U'd,12)  25+1L*(0,i)  b  2p  +  1  L  (0,1) 

which  we  obtain  from  (2,8)  by  elementary  calculations.  Prom  (2.3)  and  (2.8)  we  obtain 


(A.d. 12) 


<  const  , 


<  const ( - - - r  +  S4)  |“TT  4  conBt(b  -  1  +  6  ^ 


,  , - _= - i  m  _  *  const ;  b 

L*(0,x)  2 ( 2p  +  1)  L  <0,X)  (b  -  ’> 


„  /  D 

Ip1 1  <  const(,~ - j 

t "(0,1)  <b  -  »> 


'  -  -  -  *“==* - a*  -  _  <  const  (—■ S-  -  + 

L  <0,x)  (2p  ♦  V  L  (0 ,x)  (b  -  I)2  b  +  1  b  ~ 


<  const('g-j— -  +  6°)  , 


*p"l  m  <  const( - 

L  (0,1)  (b 


2  +  ^  )( - j  *  6*)  <  const( - — -  ♦ 

-I)2  (b  -  I)2  (b  -  I)4 


The  assertions  on  ^(x)  -  -3*( 0 )  HU  -  *> -  renain  tQ  ^  shown 

(b  +  1 )2 ( 2p(x)  +  1) 

‘P/  .  <  const [ 5(^(0)  l/i  -  x  < 

1  r  /n  d+1  b  +  1  « 

L  (0'1>  L  (0,1) 

<  con3t(~-'7)2  +  4‘ i]^  (rH-}2  4  “n8t(^ + 


|p*(x)|  -  |5*(0) 


- ^ - (i  )|  < 

2(b  +  1 ) ( 2p( x)  +  1)  (2p  +  I)2 


«  const  |  J 1 1|<(0 )  l-j-1  | - 1 - 1 

2p(x)  +  1 


by  (A.d.11).  Thus, 


‘  ~“t(^  *  <^4)2  <  *  ‘ 


'Pi1  »  <  const (- — - —  +  «4/tb/b)  , 

L  (0,1)  /b  -  1 


<pt;i  m  -  lj2*(0) 
1  L  (0.1) 


- ^4 -  h-3n_-  x)S  } 

(b  ♦  1)(2p  +  1)  2( 2p  +  I)2  L  (0,1) 


<  constl - — - -  +  (58  l/ln/K  —  ~ _ 1  £ 


APPENDIX  Bl  ESTIMATES  OP  THE  RESIDUAL 


Substitution  of  the  formal  asymptotic  approximation 

y  •  (♦  ♦  ♦  +  1^  +  Xi)^,  n  +  n  +  Xn1  +  Xn^  p  +  p  +  Xp1  +  Xp^  J  +  xj^ 

into  problem  (3.1a)  gives  the  residual 

Pxj6#b(y)  -  U2}"  +  X3}“,  -  y  (n  -  n(0)  -  x5'(0))f  +  (n(0)  -  5)^  + 

+  (n1  (0)  -  n^  )V  +  (4>*  1 0 )  -  $')n  -  X^'n^  +■  ip'n  ♦  (n^  +  n1  !  - 

-  X2^^  +  n^,  y  (p  -  5(0>  *  xp'(0))i  +  (p  -  p(0)>^  + 

e  s 

+  (p,  - 5,(0))*  +  (**  -  ?'(0))p  +  x[*,p1  +  +  (5,  +  pn )♦ 1 J  + 

♦  *2^(5,  +  PjKO.ia-1)  +  X^tX^J.O.nU"1)  +  Xn1(X"1),p(X"1)  +  Xp^X-1)) 


We  apply  Lemma  A.d.20  and  obtain 


2 

IX2^"*  „  <  const  X2( - - - —  +  S8)  , 

L  (0,1)  (b  -  1)’ 


IX2*- 1  _  <  const  X2( - - — —  +  68/ tn/b)  . 

L  (0,1)  (b  -  1)5/2 


We  now  apply  Taylor's  theorem  and  the  Lemmas  A.d.6  and  A.d.20. 


•4-  (n  -  5(0)  -  x5'(0))i«  <  y  «5  -  5(0)  -  x£'(0)i  ^  XiJ-i  + 

L  (0,1)  L  ( 0 ,4Xtq  )  l’iO,-) 

*  y  *5  “  5(0)  -  xp'(0)i  „  •  I4>i  „  < 

L  (0,1)  L  (4T0,«) 


-65- 


<  const[ At* Ip" I 


L  (0/4Xtq) 


|D<(0)|  +  ipi  „  «x3*(o))  < 

L  (0,1) 


<  const  X2(in  +  68{4n/b)2J 


Similarly  we  obtain 


I  ( n ( 0 )  —  n)4»  I 

L  {0'1>  '  l“(o,6Xt0)  t 1  l\o.-) 


Xl*,l 


L  (0,1)  L  (6TQ,») 


<  6A  T  Ip*  I 


O'**  '  -  I't'-tO)  | 

L  (0,6Atq)  1 


<  lp(0)  -  pi 

+  lp(0)  -  pi  ^  Id  I 

r  /  ft  <  \  I 

const  X3(i  in  ~  *  64/ b  (ln/b)5/2) 

<  const  X2  in  {  /wi  (g-JLj-  +  +  fi4/^  < 


+  Ip' I 


L  (0,1) 


<  const  A2  in  y  (g-  +  «8in/b)  , 


X 


.(n-^O)  -  n-,)*'  <  3XVp;1  ^  A.*l  + 

L  l0'V  L  (0,3Xx  J  l\o,»)  1  L  (0,1)  l“ 


(3V") 


<  const [ A2 in  {  /<WE  +  fiVCi)in/b  +  ( - 1—  +  S4/wi)x3/WbJ 


✓b  -  1 


<  const  A2  in  j  64(in/b)2  , 


K*-(0)  -  *')nl  ,  <  3Xxoiri  _  Xlnl  +  |*»,  ,'nt 

L  (0,1)  L  (0,3Xtq)  l  (0,»)  l"(0,1)  l”(3t  ,») 


<  const [X2 in  -1  /in/b  (g  [  —  +  «8)/in/b  +  ( - £ -  +  S8)X3  b 

(b  -  I)4  b  '  1 


in/b] 


<  const  A2  in  -J-  f-1  +  i8ln/h) 


<  const [X2 (.J- ~ — ^  ♦  «4)tn  X  (}==  +  (in/¥)5/2)  + 

+  x3[__b_  +  6*){—L-  In  {  ♦  «Vb  (in/b)5/2)J  < 

(b  -  I)2  /b  -  "1 

<  const  X2tn  {  Og-rT  +  s8/*  (ln,/b)3)  ' 


I  Xip 


•nl  <  X2lill  _  *"•  , 

1  x,1  { 0 , 1 )  L  (0,1)  L  (0.«) 


<  const  x2(^==  +  5  /in/b)  /Wb 


<  const  X^(«g  +  $  tn/b) 


,  *  *  *!,E-  *  * 


+  A I  n,  +  n.l  *  ^  i  *  a 

1  1  L  ( 4Xt  ,  1 )  L  (4T  ,-) 


f  >  ttlD 

*  const [A  i - g 


2</b  “  -  +  64/*n/b  + - b_  + 


/b  -  1 


♦  iVb  <ln/b|5/2)tn  1  <Ci  (j  ♦  S*/E  (*»^|5/J) 

„  xJ(!5zi  ,  ,vwi .  »n-j- «.  { ♦ «4^  <«»/F,5'2))(i  *”  t  * sV5: 

*■  b  /T^7 

<  const  X2*n  j  (£  +  68bUn/b) 1  V2)  . 


2-  -  .  -  ^  X  ^  ■  iD  •  I  f  ■  p  *  I  +®ni®«o  )  * 

'X  **,(»,  +  n1>'L1(0(1)  **1  1  L* (0,1)  1  L  (0,-) 


<  const  x2(^;  ♦  sV^H^1  -  +  -==  +  (tn/b,5/2)  4 

/b  -  1 


.2,  1 

<  const  X  It— — 


+  58/b  (tn/b)2) 


The  remaining  terms  are  treated  analogously,  and 


(B.1) 


,FA,«,b<y)'2.l 


<  const  A^[ - - - —  ♦  In  4-  S^lin/F) 1  ^2] 

(b  -  1> 


follows. 
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